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Abstract. It is shown here that concurrence between advection and diffusion in a drying sessile drop of a 
biological fluid can produce spatial redistribution of albumen and salt. The result gives an explanation for 
the patterns observed in the dried drops of the biological fluids. 



PACS. 47.57.-s Complex fluids and colloidal systems 
66.10.Cb Diffusion and thermal diffusion 



47.55.nb Capillary and thermocapillary flows 



1 Introduction 



The dried droplets of the biological fluid have very com- 
plex structure. In particularly, there are protein ring at the 
periphery and dendritic crystals in the central area of the 
sample (Fig.QJ. The visual appearance of a sample is used 
for diagnosing a wide range of diseases. Rather clear re- 
lationships of the 'pattern — pathological processes' type 
have been revealed 
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Fig. 1. A dried drop of the multi-component fluid consisting 
of NaCl (0.9 %), albumen (9 %), and water. This composition 
and the concentrations are typical for the blood serum of the 
healthy persons. 



The priority in the investigation of the dehydration 
self-organization phenomenon in biological fluids belongs 
to Rapis j2j. While seeming simple, the effect turned out 
to be extremely complicated and to involve a number of 
interrelated processes of different nature. 

The biological fluids are complex colloidal systems. 
Even the slightest variation in the composition of a liq- 
uid leads to total changes in the dynamics of phase tran- 
sitions during drop drying 3 . Despite the application of 
the effect to practical medical diagnostics Qj , and the dras- 
tic improvement of our understanding of pattern forma- 
tion, complex non-linear dynamics of pattern formation in 
drops of biological fluids is mostly unclear. Nevertheless, 
the performed analysis allows to conclude that the main 
effects observed in the dehydration of biological fluids are 
typical for colloidal solutions in general and can be de- 
scribed in the framework of conventional approaches 0]. 

This work is an part of the research on mathematical 
modeling of pattern formation in a biological fluid upon 
desiccation, with a special emphasis to evaporation, capil- 
lary flow , effect of diffusion '6"7I and sol-gel phase tran- 
sition. The main scientific goal of the research is determi- 
nation the relationship between the pathological processes 
occurring in an organism, the variation of the physical and 
physicochemical properties of the biological fluids caused 
by these processes, and the kind of the patterns produced 
in the drying of the sessile droplet of a biological fluid. 
The results presented in this paper allow to give the fun- 
damental explanation of spatial component redistribution 
as a consequence of diffusion and capillary flow. 



2 Model 

2.1 Theory of solute transfer 

Detailed investigation of of solute transfer in the desic- 
cated drops was performed by Deegan et al. "H*. In their 
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Fig. 2. Schematic of relevant parameters for the theory. 



theory, an outward flow in a drying drop of liquid is pro- 
duced when the contact line is pinned so that liquid that 
is removed by evaporation from the edge of the drop must 
be replenished by a flow of liquid from the interior. This 
flow is capable of transferring 100% of the solute to the 
contact line and thus accounts for the strong perimeter 
concentration of many stains. Furthermore, the theory re- 
lies only on a generic property of the substrate-solvent 
interaction, the presence of surface roughness or chemi- 
cal heterogeneities that produce contact line pinning, and 
therefore it accounts for the ubiquitous occurrence of ring- 
like stains. 

Indeed, as soon as evaporation begins, particles deposit 
onto the substrate, resulting in a strong anchoring of the 
three-phase line. In the case of biological fluid, proteins are 
adsorbed on the substrate, leading to a strong anchoring 
of the triple line. 

Recently, segregation in multi-component ceramic col- 
loids during drying of droplets was investigated [Sj. 

If the size of the solute particles is small, diffusive cur- 
rents become comparable to the advective currents [S]. 

Biological fluid widely used in the medical tests is 
blood serum. It consists of the proteins (mostly albumen) — 
about 9 %, NaCl — 0.9 %, and water. It is shown in 
this work, that concurrence between diffusive and capil- 
lary flows can produce more uniform distribution of salt 
along a drop diameter, but has not essential effect on al- 
bumen distribution. 

The conservation of fluid determines the relationship 
|H] between the vertically averaged radial flow of the fluid, 
V, the position of the air-liquid interface, h, and the rate 
of mass loss per unit surface area per unit time from the 
drop by evaporation, J: 




(1) 



where p is the density of the liquid. All quantities are 
supposed to be dependent both of distance r, and of time 
t. For additional explanation see Fig. |21 

The flow velocity can be expressed from by rewrit- 
ing it in the integral form 151 
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The conservation of solute gives 

d ^ d 

— (pch) + --^ (rphcv) = 0, 
at r or 



(3) 



where c is concentration, v = v -I- Vd is velocity including 
both advective term, v, and diffusive one, Vd, jJOl- Here 
cvd = — I?gradc, where D is the diffusion constant for 
solute in solvent. Then 
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Master equation describing spatial and temporal solute 
dynamics is 
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(5) 



One needs to know the evaporation law, the drop pro- 
file, and the vertically averaged radial flow of the fluid 
inside the drop to solve the master equation. 



2.2 Simplifying assumptions 

We will suppose that 

— air-liquid interface has a spherical cap shape induced 
by surface tension [H|; 

— drop apex height steadily decreases with constant ve- 
locity vq: /i(0, t) — ho — iigt, where ho is initial height 
of drop apex [TJj; then J = pvo/2; 

— in a multicomponent fluid the components do not in- 
fluence on each other; 

— diffusion coefficient for each component D is constant; 

— variation of the solution density p vs. concentration 
c is rather small, so we will suppose constant density 
everywhere, but the term containing dp/dc, because 
this term may be very large even for small variations 
of the density. 

Detailed description of the main assumption can be 
found in Ref. 



2.3 Evaporation models 

We could not find any experimental data about spatial 
distribution of vapor flux near the drop surface, the theo- 
retical models are rather controversial and contradictory. 

The analogy between diffusive concentration flelds and 
electrostatic potential fields (they both satisfy Laplace's 
equation) is widely used for calculation of the vapor flux 
of a sessile droplet [8113114) . In particularly, this approach 
was used to describe the vapor flux of very thin evaporat- 



ing droplets |12ll5j . In this model, the evaporation flux 
is almost uniform all over free surface except very narrow 
region near the drop edge Nevertheless, there is a 
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singularity at the contact line. The divergence at the edge 
is nonphysical and can be removed by using a smoothing 
function ,15j . 

An alternative method for determining the evaporative 
mass flux is to use heat transfer analysis ^Bl- Therefore, 
this form cannot be applied to a droplet in which the con- 
tact line is pinned because the presence of colloidal parti- 
cles affects the evaporation near the edge To mimic 
this effect the evaporation function was modified |17| . 

On the other hand, there is a region near the contact 
line with high evaporation rate |18| . 

Following Ref. JHI the simplest possible kind of va- 
por flux is assumed in the present work, i.e. a constant 
evaporation rate all over free surface. 



of albumen concentration p(c) 1000 + 300c kg/m'^ and 
iff = 0.3 1111. 

Equation ijHJl was solved numerically. Presented in Fig.|Sl 
results show that diffusive processes can prevent carrying 
out of salt on edge by capillary current. At the same time, 
diffusion has not essential influence on spatial distribu- 
tion of albumen. As a result the edge of a sample con- 
sists mainly from albumen, while its central part is com- 
posed both from salt and from albumen. This result con- 
firms qualitative conclusions obtained earlier using sim- 
plest model (6I7| . One can see that near the drops edge 
there is rather narrow region with high concentration of 
albumen. Phase transition of albumen from sol to gel have 
to occur in this region. 



2.4 Dimensionless form of the master equation 

Using above assumptions, one can rewrite master equation 
((^ in the dimensionless form 



dc _B d f dc 
dr X dx \ dx 



dc [DdL Odpdc \ kc f diy 

where B = hi^D / {vqR^) is dimensionless parameter, t = 
tvo/ho is dimensionless time, x = r/R is dimensionless 
distance, and k — ho/R. L(x, r) is dimensionless profile of 
the drop written as 



L{x,t) 



'A^ + 1 
2A 



2A ' 



where A{t) ~ k{l — r) and 

B{x, t) = ^{A^{T) + lf~{2A{T)xf ~ A^{t) - 1. 
The dimensionless velocity can be written as 



v{x, t) 



4A{t)xL{x,t) 



3 Results and discussion 

In medical diagnostics, drops of a biological liquid in vol- 
ume 10-20 ^1 are used. The drops are placed onto the 
strictly horizontal glass. Thus diameter of a droplet is 5- 
7 mm and ho fa 1 mm |T|. Hence, k — 0.33. Diffusion 
coefficients are Da — 7.7 • 10^^^ m^/sec (albumen) and 
= 1.5 • 10~^ m^/sec (NaCl). Experiments on measure- 
ment of height of a drying up drop show, that wq « 4 • 
10"'' m/sec Then dimensionless parameters are Da ~ 
0.03 (albumen) and Dg « 0.6 (NaCl). We suppose that salt 
has no effect on solution density because of very small 
concentration. Then density of solution is linear function 
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